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O ■ Abstract 

(N , 

£H \ We show that a previously proposed cosmological model based on general relativity 

with non vanishing divergence for the energy-momentum tensor is consistent with the ob- 
served values for the nucleosynthesis of helium for some values of the arbitrary parameter 
a presented in this model. Further more values of a can be accommodated if we adopt 
the Randall-Sundrum single brane model. 
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Motivated by the desire to seek a solution to the entropy problem within standard classical 
cosmology we have previously proposed a simple modification of GR; namely relaxation of the 
^ ; condition of the vanishing of the divergence of the energy-momentum tensor [1] 

O: T^;// = (1) 

This can be motivated by noting that the assumptions that lead to (JlJ are all questionable, 
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including the principle of equivalence [2]. Relaxation of (JIJ does not upset the success of 
general relativity either in or outside cosmology. It does, therefore, appear that the covariant 
conservation has not been specifically tested by observation. The construction of a theory 
in which (JIJ does not necessarily hold would thus provide an opportunity for testing this 
condition. In fact the relaxation of this condition, introduces a single arbitrary constant, a , 
where < a < 1. The field equations in this model become [1]: 

- ^Rg^ = -KT^ (2) 

Where: 

7= o 7T ' K = • ( 3 ) 

S — 2a a 

Imposing the Bianchi identity one obtains: 

T"";/i = -i(l-a)I> (4) 

Standard GR has a = 1. 

The vacuum field equations remain the same as in standard GR. Consequently the crucial 
tests of GR and the important analytic features (such as the existence of singularities and 
black holes) are maintained. The Field equations (0) are also the same as the standard field 
equations for systems with traceless energy-momentum tensors, except that G is replaced By 



G/a. This shift in the effective value of G to G/ot in the early universe is expected to change 
the prediction of the primordial nuleosynthesis. In this paper we use a simple, approximate and 
semi analytical method to check the consistency of this model with the observed values for the 
nucleosynthesis of helium. We show that at least for some values of a the model is consistent 
with the observed values. This problem was investigated in ref. [3] using a different approach. 
The modified filed equations, in the radiation dominated era in this model become: 
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This represent an increase in density as related to standard model. 

The primordial production of 4 He is controlled by a competition between the weak inter- 
action rates and the expansion rate of the universe. As long as the weak interaction rates are 
faster than the expansion rate, the neutron-to-proton (~J ratio tracks its equilibrium value. 
Eventually as the universe expands and cools, the expansion rate comes to dominate and ^ 
essentially freezes out at the so called freeze out temperature. However the nucleosynthesis 
chain which begins with the formation of deuterium through the process p + n D + Y is 
delayed past the point where the temperature has fallen below the deuterium binding energy 
Eb since there are many photons in the exponential tail of the photon energy distribution with 
energies E > Eb despite the fact that the average photon energy are less than Eb ■ (see for 
example [4]). 

The increase in density is expected to modify the expansion of the universe during nucle- 
osynthesis by modifying both the freeze out temperature and the time available for the decay 
of neutron. These two factors will lead to an increase in the nucleosynthesis of Helium. 

Equation @ can be written as : 

R\ 2 8nG 8nGl-a 
RJ 3 3 a 

The perturbation in density presented in Equation. (jUJ) can be taken as a modification of 
the gravitational constant 5G, where 

1 - a 

5G = - — -G (7) 

a 

Use of the conservation equation for radiation with equation (JEJ) leads to: 

t s 

Where t m is the modified time from freeze out to start of nucleosynthesis and t s is the time 
calculated from Standard Model without the extra energy. 

The calculation of primordial element abundances is a highly nonlinear problem with many 
coupled nuclear reactions, and requires a numerical analysis. There are many different numerical 
codes for doing this calculation starting with Wagoner [5]. They mainly differ in the different 
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factors they include in their calculations and particularly the different values they use for the 
neutron half life [6, 7]. Here we only present a simplified and approximate method. In this 
approximation, the primordial helium abundance Y p is given by: 



^=(jv^r) y exp[ - A(i "'- iF)1 



Where F represent freeze out and we ignored tp because it is of two order less than t m 
Equation (JOJ) can be written as 

Y = ( J exp( — -pJ-) where x = (10) 

\l + xJ F Va N p 

And thus 

„ — N n 
X N p 

mined by equating the Hubble constant H to the weak reaction rate rj for n <-» p conversions. 
Now, 

H oc G 1/2 Tj w and rj cc T / 5 io 



exp (jr^) , where T/ 10 is the freeze out temperature MT n = j^r^j and is deter- 



So that: Tf io = G 1 ' 2 



- — 

(x\nx)e ^ 5G 

■ 6Y = - 3(1 + *)* ~~G (12) 
which from equation (J7J) becomes: 

(xlnx)e 1-q 

5F = - 3(1+*)* ( ~ } (13) 

Putting x = 0.14 (From Y = 0.246 [6]), A = 78 x 10~ 5 , t s = 120s, and if we take a value for 
a = 0.9, equation (JT3*|) gives: 

SY = 0.0071 

which is quite within the margin of difference between measured and observed values. This 
value, as expected will increase, for smaller values of a . 



3 



Discussion 

Big bang nucleosynthesis (BBN) is one of the most sensitive available probes of physics 
beyond the standard model. The 4 He abundance in particular, has often been used as a 
sensitive probe of new physics [8 — 12]. This is due to the fact that nearly all available neutrons 
at the time of BBN end up in 4 He and the neutron-to-proton ratio is very sensitive to the 
competition between the weak interaction rate and the expansion rate. Observations of light 
element abundances have improved dramatically over the past few years. The recent all-sky, 
high-precision measurements of microwave background anisotropics by WMAP has opened the 
possibility for new precision analysis of BBN. Using BBN prediction gives a powerful constraint 
over the various cosmological parameters. The possibility of the physical 'constants' taking 
different values at different times in the universe history has recently received much attention 
with the apparent observation that fine structure constant had a different value in the distant 
past [13]. Variation of the gravitational constant, that was originally started by Dirac [14] in 
his so called " Large Number" theory is now constrained by the observations of the light element 
abundances [15-18]. 

Copi et al [16] uses the recent measurements of the primordial deuterium abundances (D/H) 
in conjunction with the WMAP determination of rj to set a Limit for G/G , where G is the value 
for the gravitational constant during nucleosynthesis and Go its value at present. They found 
that G/Gq is constrained in the range 1.21 to 0.85 at the 68.3% confidence level and between 
1.43 and 0.71 at the 95% confidence level Assuming a simple power dependence G ~ t~ x , x was 
constrained to the range— 0.004 < x < 0.005 at the 68.3% confidence level, —0.009 < x < 0.01 
at the 95% confidence level. In our model taking a = 0.9 gives G/Gq = 1.11, which is within 
the range obtained by Copi et al .Further Cyburt et al [17] using 4 He abundance set a more 
restrictive limit for AG/G of 13%. 

Finally we mention that if one adopted the recent ideas of brane model, where the universe 
is supposed to be embedded in a higher dimensional bulk, particularly the Randal Sunundrun 
type II model [19], one could accommodate more smaller values for a .The predicted expansion 
rate in the early universe can be significantly modified, particularly with the presence of a dark 
radiation term [20,21], which may take negative values and can contribute as high as 27% of 
the background photon energy density [22]. 
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